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Aims. We study the LMC interstellar medium in the field of the nebula N157B, which contains a supernova remnant, an OB associ- 
ation, ionized gas, and high-density dusty filaments in close proximity. We investigate the relative importance of shock excitation by 
the SNR and photo-ionization by the OB stars, as well as possible interactions between the supernova remnant and its environment. 
Methods. We apply multiwavelength mapping and photometry, along with spatially resolved infrared spectroscopy, to identifying the 
nature of the ISM using new infrared data from the Spitzer space observatory and X-ray, optical, and radio data from the literature. 
Results. The N157B SNR has no infrared counterpart. Infrared emission from the region is dominated by the compact blister-type 
HII region associated with 2MASS J05375027-691 1071 and excited by an 08-09 star. This object is part of an extended infrared 
emission region that is associated with a molecular cloud. We find only weak emission from the shock-indicator [Fell], and both the 
excitation and the heating of the extended cloud are dominated by photo-ionization by the early O stars of LH 99. 
Conclusions. Any possible impact by the expanding SNR does not now affect the extended cloud of molecules and dust, despite the 
apparent overlap of SNR X-ray emission with infrared and Ha emission from the cloud. This implies that the supernova progenitor 
cannot have been more massive than about 25 M Q . 
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1. Introduction 



The Ha-emitting region N157 B ( Henize| 1956 1 in the Large 
Magellanic Cloud is on the southwestern fringe (projected lin- 
ear distance 90 pc from the center) of the major star-forming 
complex 30 Doradus. In a limited area of (projected) diame- 
ter 65 pc, the field contains a supernova remnant, an OB asso- 
ciation, HII emission regions, bubbles of low-density hot gas, 
as well as the neutral material in dense clouds (Laz endic et al.| 
(2000) and references therein). The structure, dynamics, and en- 
ergy balance of this region are poorly understood. There is even 
uncertainty about the nature of the emission observed at different 
wavelengths. 

The OB association LH 99 ( |Lucke & Hcdg1p9"70"] l is seen 
on the sky in the same direction as N157B, and is often consid- 
ered to be associated with it because it has the same foreground 
extinction (Wang & Gotthelf 1998|l. LH 99 c ontains a large num- 
ber of (early) O stars ( |Schild & Testor|1992[ ), powerful ionizing 



SNR emission. However, various observations strongly support 
the latter interpretation: bright X-ray emission coincides with the 
radio source ( Long & Helfan d*|i979), th e ratio of radio-to-optical 
intensity is quite high ( Dickel et al.|1994) , spectral line images 
( |Danz iger et a l.|1981| l reveal filamentary structures with a high 
line ratio[SH]/Ha > 0.7, the lines of [OI], [Fell], and [Felll] are 
strong, whereas the Hel /14686 A line is weak. The SNR has a 
peculiar one-sided morphology, does not have a well -determined 
outer boundary, and exhibits bright radio/X-ray core surrounded 



by an extended envelope of linear size of 30x20 pc ( Chu et al. 
[T9921 |Dickel et al.|[l994l |Lazendic et alj|2"000] i, or even larger 



( |Townsley et al. 2006| l making the remnant unusually large in 



sources whose strong stellar winds are expected to produce bub 
bles of low-density hot gas ( Townsley et al.|2006 1 and to illumi- 
nate the dusty filaments and clouds such as the nearby molecular 
cloud 30Dor-22 ( [Johansson et al.|1998) . 

Much of the nebula N157 B can be identified with the Crab- 
like supernova remnant (SNR) B0538-691. Its embedded X- 
ray pulsar PSR J0537-6910 suggests an age of about 5000 yr 
( |Marshall et al.|1998| |Wang & Gotthelf||1998|>. The N157B ra- 
dio counterpart MC 69 ( |Le Marne||1968| |McGee et al.111972 
has a spectral index a = -0.19, S v oc V* ( |Mills et al.||l978 
|Lazendic et al.| 2000] >. Such a flat spectrum is characteristic of 
thermal HII region emission, as well as of nonthermal Crab-type 



both radio and X-ray emission ( Lazendic et al.||2000| [Wang & 
|Gotthelf|1998) . 

In such a complex environment, shocks from SNRs and 
photons from luminous stars compete as heating agents for the 
ISM. The conditions in and around N157B have therefore been 
a frequent subject of study. Several authors, e.g. |Chu et aL] 
( f!992> ; |Townsley et~al] ( |2(5o6> ; |Chen et al.| have con- 

cluded that the expanding SNR is presently colliding with the 
molecular cloud to its south. The impact of a shock front on 
dense cloud material should leave a signature detectable at mid- 
infrared (MIR) wavelengths ( |Ohva et aL]|1999| |Reach & Rho 



2000). Confirmation of the proposed collision is of great interest 



Send offprint requests to: E. R. Micelotta 



as it would identify an excellent source for studying the physics 
of ongoing dense ISM processing by shocks. 

The use of infrared observations is particularly well-suited 
to trace physical conditions in dusty environments because of 
their insensitivity to extinction. In the following, we combine 
our Spitzer observations with literature data at a variety of 
wavelengths to study the role of the different physical compo- 




Fig. 1. IRAC and MIPS images of the N157B region, calibrated in MJy/sr. The diamond marks the peak of the supernova remnant 
X-ray emission (Sasaki et al. 2000| l and the circle marks the 2MASS compact source J05375027-69 11071. The coordinates are (RA, 
dec) J2000. 
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Fig. 2. Left panel: IRAC 8.0 pm map overlaid with Chandra ACIS-S X-ray intensity contours of N157B at 4, 5, 7, 8, 12, 16, 20, 
25, 30, 35, 50, 100, 500, 30000 x 10~ 2 counts sec" 1 arcsec" 2 (courtesy D. Wang). The rectangles mark the IRS slit locations: SL = 
east-west black, LL2 = north-south black, SH and LH are the two small white boxes in the center of the field. The 03 stars ST 1-62, 
1-71, 1-78 in LH 99 ( Schik f~& Testor|1992 1) with projected positions inside the N157B X-ray contours are marked with an "X" 
( |Schild & Testor|1992| >. Right panel: MCELS Ha image of N157B ( |Smith & MCELS Team|T998] ), again with the 03 stars marked. 
Note the correspondence between the 8 pm IR emission and Ha extinction. The compact object J05375027-691 1071 is marked by 
a circle which also identifies the area over which we determined the Ha/[S II] ratio (see Sect. 4.1). In both panels the coordinates 
are (RA, dec) J2000. 



nents (SNR, HII region, OB association, dust clouds) and heat- 
ing processes (shock-excitation, photo-ionization) in shaping the 
N157B region. Throughout, we adopt a distance of 50 kpc for the 
LMC (I Westerlund|1997[ ), so that 10 " correspond to 2.4 pc. 



2. Observations and data processing 

We have observed N157B with the Infrared Array Camera 



Table 1. IRS observational setup. The IRS modules are SL: 
Short-Low, LL2: Long-Low 2, SH: Short-High, LH: Long-High. 



(IRAC ) (|Fazio et al.||2004) > and with the Infrared Spectrograph 
(IRS) ( Houck et al.|2004[ )~on board the Spitzer Space Telescope 



(Werner et al. 2004| l, as part of the IRS guaranteed-time 
program (PID 63, PI J. R. Houck). The images from the 
Multiband Imaging Photometer for Spitzer (MIPS) have been 
retrieved from the archive (MIPS 24 pm: program 3680, PI K.J. 
Borkowski, MIPS 70 pm: program 20203, PI M. Meixner). 

The IRAC data (Fig.[T]i were taken on 7 December 2003 and 
consist of mosaic images in four channels at 3.6, 4.5, 5.8 and 8.0 
pm. The Basic Calibrated Data (BCD) products from the Spitzer 
Science Center (SSC) pipeline were used to construct the mosaic 
images. 

The IRS spectra were taken on 26 May 2005 using the stan- 
dard IRS "Staring Mode" Astronomical Observing Template 
(AOT). The observational setup is reported in Table [T] The IRS 
slit coverage is shown in Fig. [2]- left panel, outlined with rectan- 
gles: SL goes from East to West, LL2 from North to South and 
SH and LH are the two small rectangles in the center of the field. 

Each cycle yielded two exposures at different nod positions 
along the slit. The data have been pre-processed by the SSC data 
reduction pipeline version 12.0.2 (Spitzer Observer's Manual, 



Cycles 


Duration" 


IRS 




Resolving Power 


number 


(sec) 


module 


(pm) 




2 


14 


SL 


5.2 - 14.5 


64 - 128 


2 


14 


LL2 


14.0-21.3 


64 - 128 


3 


30 


SH 


9.9- 19.6 


600 


3 


60 


LH 


18.7 - 37.2 


600 



(a): Total integration time: SL + LL2 = 1 12 sec, SH + LH = 540 sec. 



chapter The two-dimensional BCD constituted the basis for 
further processing. First we corrected the BCD frames for bad 
and rogue pixels with the IRSCLEAN^tool. IRSCLEAN identi- 
fies bad and rogue pixels and replaces them with the average of 
good nearest neighbors. Rogue pixels are single-detector pixels 
that show time-variable and abnormally high flux values. Next, 
we combined the frames from the same nod position using the 
mean where two frames were available, and using the median in 
case of three available frames. For each integration at the same 
position on the target, an off-source sky exposure is provided, 
so we computed the mean/median sky (depending again on the 
number of available frames) and subtracted it from the corre- 
sponding nod position frame. 

Further processing was done using the SMART package ver- 
sion 5.5 (Higdon et al. 2004 1 - a suite of IDL software tools 



http://ssc.spitzer.caltech.edu/documents/SOM 
available at http://ssc.spitzer.caltech.edu/archanaly/contributed/ 
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Fig. 3. Composite image of the N157B region: red — Spitzer 
IRAC 8.0 fxm, green = MCELS Ha, blue = Chandra *ACIS* 
0.9 - 2.3 keV. The blue area in the center shows X-rays from the 
northend parts of the remnant filling a prominent hole in the 8.0 
pm. The contours of the optical and infrared emission (see Fig. 
[2]) are recognisable in the cyan and magenta areas respectively . 
This picture is part of a large image of 30 Dor from Townsley et 
al. (2006). (Reproduced by permission of the AAS.) 



Table 2. Infrared flux densities. 



Origin 


A 


Bandwidth 


Flux Density (ly) 




(pm) 


(pm) 


Compact 
Object 


Extended 
Cloud 


2MASS 


1.24 


0.16 


0.0013 ± 0.0001 






1.66 


0.25 


0.0022 ± 0.0002 






2.16 


0.26 


0.017 ± 0.006 




IRAC 


3.6 


0.75 


0.13 ±0.02 


1.7 ±0.3 




4.5 


1.01 


0.33 ± 0.05 


3.2 ±0.5 




5.8 


1.42 


0.66 ±0.1 


7.2 ± 1.1 




8.0 


2.93 


0.98 ± 0.2 


19 ±3.0 


MIPS 


24 


4.7 


1.7 ±0.3 


38 ± 6.0 




70 


19 


3.6 ± 0.5 


139 ± 21. 




160 


35 


2 ±0.3 


105 ± 16. 


IRAS 


12 


7.0 




3 ±0.5 




25 


11.2 




22 ± 3.0 




60 


32.5 




124 ± 19. 




100 


31.5 




312 ±47. 



developed for spectral extraction and spectral analysis of IRS 
data. The high-resolution spectra were extracted using full slit 
extraction. The low-resolution spectra were instead extracted us- 
ing SMART'S interactive column extraction. The calibration is 
based on observations of standard stars ( |Decin et al.|2 004). The 
ends of each orders, where the noise increases significantly, were 
manually clipped. To obtain a good match between the low and 
high resolution modules, the SL + LL2 spectra were scaled up 
by 10%. Finally the resulting spectra from the two nod positions 
were averaged to obtain the final spectrum shown in Fig. [4] 



PAHs 

6.2 7.68.6 11.2 12.7 




20 25 

Wavelength [^m] 



Fig. 4. Combined low- and high-resolution IRS spectrum of the 
2MASS source J05375027-691 1071. 



3. Results and analysis 

3.1. IRAC and MIPS images 

The IRAC and MIPS images (Fig. [[} of the N157B region are 
dominated by emission from the dust associated with the molec- 
ular cloud south of the X-ray emission region (peak marked by 
a diamond). Relatively faint in the 3.6 pm and 4.5 pm images, 
the irregularly shaped diffuse emission is quite bright at 5.8 pm, 
8.0 pm, and 24 pm where it appears more extended with a quasi- 
circular shape. In the 70 pm image, it no longer stands out clearly 
as it suffers from confusion with the very extended low-surface 
brightness emission characteristic of the whole 30 Doradus re- 
gion. The cloud is listed as object no. 1448 in the IRAS catalogue 
of LMC sources by Schwering] ( [19 89 ) and shows up in the rela- 
tively low-resolution IRAS maps as an extension of the main 30 
Doradus IR source. Fig.[2](left) and Fig.|3]show that the contour 
of the 8 pm emission (sensitive to PAHs) from the cloud follows 
reasonably well the outline of extinction in the Ha image (Fig. [2] 
- right). Much of the infrared-emitting dust must therefore be 
either in front of the ionized gas, or embedded within it. From 
the MIPS 24 pm, the size of the cloud is roughly 2' (28.8 pc). 
Its dimensions are thus similar to those of the SNR but the two 
are offset from one another, as shown by the superposition of 
the X-ray and 8 pm maps in Fig. [2] left-hand panel. Although 
there is considerable overlap between the X-ray and Ha emis- 
sion regions (Fig. [2] right-hand panel), there is no trace of an IR 
counterpart to the X-ray SNR emission: the two occur almost 
side-by-side. 

The bright and compact object at RA (J2000) = 
5h37m50.28s, DEC (J2000) = -69°11'07.1" is located within 
the confines of the infrared cloud. In the IRAC images, it has a 
diameter of ~3 pc and it is the brightest source in the field in all 
IRAC bands. Centimeter-wavelength radio maps by DickeletaL] 
([1994} and |Lazendic et al^[ ( |200 0) reveal weak radio emission at 
its position. Although there is ISO SWS and LWS spectroscopy 
of the extended cloud just described ( Vermeij & van der Hulst 



2002), the ISO apertures did not include this bright object. We 
have chosen this object, which is identical to the 2MASS near- 
IR source J05375027-6911071, as the reference position for the 
IRS slits (Fig.[2]- left panel). 

We have collected infrared flux densities for both the com- 
pact object and the whole dust cloud shown in Fig. [T] by inte- 
grating the emission over circles with radii 6" and 57" centered 
on RA(5h37m50.3s), DEC(-69°11'07") and RA(5h37m45.3s), 
DEC(-69°10'47"), respectively. The results are listed in Table[5] 
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Table 3. Spitzer fine-structure lines observed in J05375027- 
6911071. 



LL2 emission lines 



Line ID 


"-rest 


EP" 


Flux* 


EW" 




\pm\ 


[eV] 


[10~ 20 Won" 2 ] 


W 


Clear detections 


[SIV] 


10.51 


34.8 


2.35 + 


0.26 


10 


[Ne II] 


12.81 


21.6 


17.2 : 


t4.0 


80 


[Ne III] 


15.56 


41.0 


8.4 ± 


0.26 


40 


[S III] 


18.71 


23.3 


13.3 : 


t0.3 


60 


[Fe II] 


25.99 


7.9 


3.32 + 


0.93 


20 


[S III] 


33.48 


23.3 


55.6: 


t 1.4 


210 


[Si II] 


34.82 


8.15 


20.7 : 


fc0.8 


70 


[Ne III] 


36.01 


41.0 


1.64 + 


0.52 


5 


Marginal detections 


[Arm] 


8.99 


27.6 


1.20 + 


0.27 


4.27 


H 2 S(l) 


17.03 




0.8 + 


0.37 


6.54 


[Fe II] 


17.93 


7.9 


0.17 + 


0.19 


1.70 


[Ar III] 


21.83 


27.6 


0.74 + 


0.20 


4.00 


[Fe III] 


22.92 


16.2 


1.47 + 


0.40 


6.53 



Note. - Emission line properties obtained from single Gaussian fits in 
the IRS spectrum (Fig.H). 

(a) : EP = Excitation Potential, EW = Equivalent Width (observed). 

(b) : Quoted uncertainties are the errors from the line fit and do not 
include calibration uncertainties. 



Table 4. PAH dust emission features measured in J05375027- 
6911071. 



10.0 



SOUTH 

Continuum at 15 fj.m 



NORTH 




-50 50 

Distance from Source 1 [arcsec] 

Fig. 5. Spatial variation of continuum and emission line intensi- 
ties along the low-resolution slit LL2, calculated with the pro- 
gram PAHFIT. The vertical bar on the left mark indicates the 
representative error on fitted intensities. Spatial distances along 
the slit are measured from the position of J05375027-691 1071 
(0 arcsec). The horizontal bar indicates the size of the five-pixel 
wide extraction window for J05375027-691 1071; all other data 
points refer to three-pixel wide extraction windows. 



cent 


Flux* 


EW 


Relative'' 


[pm] 


[10- 20 Wcm" 2 ] 


[nm] 


strength 


6.2 


36.06 + 0.58 


34.1 


1.00 


7.7 


53.00 ± 0.36 


51.5 


1.47 


8.6 


41.42 ±0.54 


40.3 


1.15 


11.2 


32.77 ± 0.64 


36.7 


0.90 


12.7 


27.45 ± 0.36 


26.7 


0.76 



(a) : Central wavelength. 

(b) : The quoted uncertainties are the errors from the line fit and do not 
include the calibration uncertainties. 

(c) : EW = Equivalent Width (observed). 

(d) : Strength relative to 6.2 fim feature. 



which for convenience also yields the relevant flux densities 
taken from the 2MASS on-line data archive, and the IRAS 
database published by Sch wering] ([1989). Although all flux- 
densities in Table[2]have very small formal errors, a major un- 
certainty (easily a factor of two) arises in the separation of the 
source from its surroundings. As noted by Schwering, this is 
clearly true for the IRAS flux densities, which are very hard 
to separate from the overwhelming emission of the 30 Doradus 
complex, but a glance at Fig.[TJ shows that this problem is not 
limited to the IRAS data, but also extends to the Spitzer map- 
ping of this complex region. 

3.2. IRS Spectroscopy of J05375027-691 1071 

The combined high- and low-resolution spectrum of the 2MASS 
source J05375027-691 1071 covers the wavelength range from 5 
jjm to 38 jum, and is shown in Fig. [4] The (short-wavelength) 



[Ne II] 1 2.81 fj.m 




-40 -20 20 40 -40 -20 O 20 40 

Distance from Source 1 [arcsec] 



Fig. 6. Spatial variation of emission line intensities along the 
low-resolution slit SL. Otherwise as Fig. [5] Towards J05375027- 
691 1071 the [S IV] line is undetectable because of strong silicate 
absorption. 



low-resolution part has been extracted from the single slit seg- 
ment centered on 2mass-j05375027-691 1071 (see below), while 
the high resolution part results from integration over the full 
length covered by the SH and LH slits (see Fig. [2] - left-hand 
panel). This area includes J05375027-691 1071 but also some of 
its fainter surroundings. As a result, the spectrum in Fig. ^ac- 
curately reflects the emission from J05375027-691 1071 short- 
wards of ~20 fim. However, at wavelengths beyond 20 //m the 



6 



E. R. Micelotta et al.: Spitzer observations of LMC-N157B 



SL PAH features 




, . PAH 6.22 fjm 




0.1 r 



PAH 7.7 



0.1 - 



W 




PAH 8.6 firn 




PAH 11.2 /im 




PAH 12.7 jLtrn 



[Ne ll]/[S IV] 




[Ne ll]/[S IV] 



20 40 -40 -20 20 

Distance from Source 1 foresee] 



Fig. 7. Spatial variation of PAH feature intensities (top panels) 
compared to the [Ne II]/[S IV] ratio (same plot reproduced in 
both bottom panels for clarity) along the low-resolution slit SL. 
The vertical bar on the left mark indicates the representative er- 
ror on fitted intensities in the region between -35 and -5 arcsec. 
The oscillations are due to instabilities in the fit caused by weak 
signals. Otherwise as Fig. [5] 

spectrum includes a non-negligible contribution from the sur- 
rounding diffuse emission. 

Different ISM processes have left their mark on the spec- 
trum: emission from PAHs, absorption by silicates around 10 
/im and in a broad band between 15 fj.m and 22 /urn, various fine- 
structure emission lines and steeply increasing continuum emis- 
sion from hot dust. The dust PAH features have been modelled 
with the program PAHFIT^ (Smith et al. 2007) and the results 
are presented in Table [4] PAHFIT is an IDL tool for decompos- 
ing Spitzer IRS spectra, with special attention to PAH emission 
features, and it is primarily designed for use with Spitzer low- 
resolution IRS spectra. The program is based on a model consist- 
ing of starlight, thermal dust continuum with temperature from 
35 to 300 K, resolved dust features and feature blends, prominent 
emission lines and dust extinction dominated by the silicate ab- 
sorption bands at 9.7 jjm and 15-22 //m. PAHFIT uses Gaussian 
profiles to recover the full strength of emission lines and Drude 



Available at http://turtle.as.arizona.edu/jdsmith/pahfit.php 



profiles for dust emission features and blends. The fine-structure 
lines have been measured by single Gaussian fits within SMART, 
and the results are listed in Table 3] All these lines, except the 
iron lines, are typically associated with ionized gas in photon- 
dominated regions (PDRs - e.g. Martin-Hernan dez et al.|[2002 
|Tielens|2005D . 

3.3. IRS spectroscopy of the extended dust cloud 

As shown in Fig. [2] the SH and LH slits cover only a very limited 
part of the cloud, mainly the compact 2MASS object. The two 
much longer low-resolution slits LL2 and SL sample a much 
larger part of cloud, even extending into the cloud surroundings. 
The slits are almost perpendicular to one another and overlap at 
the position of the compact object 2MASS J05375027-691 1071. 

In order to study the spatial variation of the spectral features 
along the two slits, we sub-divided the region covered by each 
slit into 33 segments and extracted a spectrum from each seg- 
ment. We have chosen overlapping extraction windows of three 
pixels width, moving them along the slit in single-pixel steps. 
Thus, adjacent extractions are not independent of one another 
as they overlap by two pixels. The corresponding spectrum is 
therefore effectively a boxcar-smoothed spectrum. These partic- 
ular choices resulted from tests performed on the data to find the 
minimum extraction requirements needed to obtain spectra free 
of sampling artifacts. For the same reason, the region containing 
the bright object J05375027-691 1071 was treated slightly differ- 
ently. Here, we applied an extraction window five pixels wide, 
not overlapping with the adjacent pixels. As a consequence, the 
size of the extraction window for J05375027-691 1071 exceeds 
the dimensions of the intersection region common to both the 
SL and LL2 slits; it corresponds to a 7 pc x 2.7 pc rectangle. We 
used PAHFIT with its default set of lines and continuum features 
to determine the intensities of the spectral features at every slit 
position. We have plotted the intensities of fine-structure emis- 
sion lines and dust features thus extracted as a function of the 
position along the slit, expressed in terms of the distance (pos- 
itive and negative) from J05375027-691 1071 in Figs. [5] § and 

m 

The SL spectra show PAH emission in the 6.2, 7.7, 8.6, 11.2 
and 12.7 /jm bands. In determining the intensity of the latter, 
we first removed the contribution from the [Nell] 12.8 yum line. 
We note that the PAH emission from the N157B region is re- 
markably weak compared to that of other sources in 30 Doradus 
flBrandl|2008D . 

4. Discussion 

4.1. The nature of the compact object J05375027-691 1071 

The flux densities in Table[2] are not accurate enough to war- 
rant detailed SED fitting of J05375027-691 1071. However, by 
comparing various modified blackbody fits with straightforward 
integration we find an integrated flux of 5 - 1 1 x 10~ 13 W m -2 
which implies a luminosity L = 3.8 - 8.6 x 10 4 L . This lumi- 
nosity corresponds to that of a Bl-09 star if all stellar photons 
are converted to IR emission; it thus places a lower limit on the 
spectral type of the exciting star(s). Much of the emission must 
arise from hot dust, with temperatures between 180 and 350 K. 
The exciting star must therefore be close to the dense neutral 
material, and we would expect the interface between star and 
neutral material to consist of very dense ionized gas. 

The Spitzer spectrum of J05375027-691 1071 (Fig. |4j con- 
tains several diagnostic fine-structure lines. The [SIII] lines at 
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18.7 jum and 33.5 //m arise from different levels with the same 
excitation energy and their ratio provides a measure of the elec- 
tron density (e.g. |Tielens|2005) . We find a [SIII] 18.7 /mi/ [SIII] 
33. 5/mi ratio of 0.24. Collisional excitation models (Alexander 
et al. 1 1 999] > place this ratio in the low-density limit, and indicate 
the presence of an ionized gas with n e « 100 cm 3 . Such a den- 
sity is not uncommon for a parsec-sized HII region (cf. Fig. 2 
in Habing & Israel 1979 1. It is also quite consistent with the 



weak radio emission (about 40 mJy at ,L13.5-13cm) in the maps 
published by |Lazendic et al.| ( |2000| i and |Dickel et aL] ( |1994 1. 
Assuming all radio emission in this direction to be thermal and 
optically thin free-free emission, we calculate an r.m.s. electron 
density < n\ > 5 = 100 - 250 cm . It is unlikely that the ion- 
ized gas and the hot dust have different sources of excitation. We 
must conclude that the bulk of the infrared line and radio con- 
tinuum emission arises from gas extended over a volume much 
larger than occupied by the dense ionized gas surmised in the 
previous paragraph. As we have not resolved structure on scales 
less than a few parsecs, this is quite feasible. 

The ratios of lines of the same species but arising from differ- 
ent ionization states reflect the degree of ionization and the hard- 
ness of the stellar radiation field. We have measured such pairs 
of neon and sulphur lines and find ratios [SIV]/[SIII] =0.18 and 
[NeIII]/[NeII] = 0.49. In a sample of HII regions with different 
metallicities in the Milky Way, the LMC, and the SMC these two 



ratios are tightly correlated (cf. Fig. 1 by Martin-Hernandez et al. 
|2002] >. Our result fits this correlation very well but the individual 
ratios are lower than those in the (very) bright LMC HII regions 
which have [SIV]/[SHI] = 0.6-1.0, and [Neffl]/[NeII] = 1.4-6.3. 
Not surprisingly, the gas in J05375027-691 1071 thus has a lower 
degree of ionization as is expected from excitation by stars less 
hot than the ionizing stars in the bright LMC HII regions. The 
photoionization models byJSchaerer & de Koter| ( |1997| l (CoStar) 
and by Pauldrach et al. (2001 ) show our line ratios to result from 
stellar radiation fields with T e ff = 33000 - 40000 K correspond- 
ing to spectral type 05-09 (cf. |Martins et al.|2002| . Under the 
same assumption of optically thin free-free radio emission, we 



calculate from the maps by |Dickel et al.| ([1994 ) and Lazendic 
et al.| ( p000) l a minimum required Lyman-continuum photon flux 
Nl = 48.96 s _1 . This corresponds to the output of an 07.5-09 
star( |Vacca et al.|1996| ). 

Bright [Fell] 17.9 /mi, [Fell] 26.0 /mi and [Sill] 34.8 /mi 
lines trace the return of iron and silicon to the gas phase follow- 
ing the destruction of dust grains by shocks ( |Qliva et al.|1999 



Reach & Rho 2000), but a strong [Sill] line is also frequently 
detected in PDRs ( [Peeters et al.||20T)2l >. In J05375027-69 11071 
[Sill] is one of the strongest detected lines, whereas the [Fell] 
lines are weak with respect to the others. The [SI] 25 /mi line, 
also expected in shocked gas ( Tielens|2 005), is missing but sul- 
phur is detected as [SIII] and [SIV], a situation characteristic for 
photo-ionized gas. 

The intensity of the optical lines [SII] AAblll, 6731 A and 
Ha /16563 A provides a reliable tool to discriminate between 
shock-excited and photo-ionized plasma ( |Long et al.|1990|l) and 
were, in fact, used to establish the nature of N157B (Danziger 
et al. 1981]). In SNRs, [SII] emission is usually stronger than 



in photoionized HII regions, where sulphur is mostly doubly 
ionized. The value of the [S II]/Ha ratio separating SNRs and 



1980 Fesen et al. 



HII regions is about 0.4 (D'Odorico et al. 
1 1985 Long et al. 1990| ). No optical spectra are available for 
J05375027-6911071, but we may estimate the [SII]/Ha ratio 
from the MCELS line emission maps. The lines are close enough 
to assume that they suffer approximately the same extinction. 



Extended Cloud 



IRAC i — i 
MIPS (-*• 
IRAS >~» 
Fit MBB 



T 4 = 230 K 
T 5 = 29 K 



Wavelength [ttm] 



Fig. 8. Spectral energy distribution of infrared emission from the 
extended dust cloud. The data points are fitted with two modified 
black-body with T = 230 and T = 29 K. The horizontal bars 
indicate the width of the IR AC (|Fazio et al.[2 004 ), MIPS ( |Rieke| 
|etal.|2 004 ) and IRAS bands (Schwerin g|1989| ). The 12 /mi IRAS 
point has been removed because of the strong 10 /mi silicate 
absorption of the compact source, which makes it impossible to 
separate out the extended cloud contribution. 



We find a ratio [SII]/Har = 0.2, which is typical for H II regions 



(Long et al. 



1990). Note that in most of N157B, this ratio ex- 



ceeds ~ 0.7 ( |Danzigeret aL|1 981) 



Finally, we note that the 9.7 /mi silicate absorption feature is 
very strong towards J05375027-691 1071, with an optical depth 
Tq v « 1.1. This implies the presence of a large column of neu- 
tral material in front of the ionized line emission region. 

Thus, the 2MASS source J05375027-691 1071 contains a 
purely photo-ionized gas of moderate density. The high column- 
density and elevated temperature of the dust suggests a blister- 
type geometry (Isra eli 1978 1 seen from the back, and a source of 
excitation located close to the interface of the HII region with 
the obscuring dense neutral material. In that case, no more than 
about half of the ionizing photons may escape, and we conclude 
from this and the local radiation field hardness that the excitation 
of J05375027-691 1071 is caused by an obscured but no longer 
embedded 08 or 09 star. 

4.2. The northeast edge of the dust cloud 

The intensity variation of the spectral features along the rather 
short SL slit is shown in Fig. [6] and [7] The SL slit runs from 
northeast to southwest. It is potentially of interest as it cuts 
across the X-ray SNR/IR dust cloud boundary (Fig. [2] left-hand 
panel). Apart from J05375027-691 1071, most of the emission in 
all lines comes from the dust cloud, outside the X-ray contours. 
The [S IV] intensities peak around +15 arcsec (3.5 pc) southwest 
from J05375027-691 1071. Towards the object itself, because of 
the low resolution of the SL slit we could not separate the [SIV] 
emission from the silicate absorption at 9.7 /mi, contrary to the 
high resolution spectrum where the [SIV] is well detected. Both 
[Ar III] and [Ne II] have intensity distributions strongly peaking 
on J05375027-6911071, and that of [Arlll] is almost symmetri- 
cal. The [Nell] distribution has an asymmetry similar to but less 
outspoken than that of [SIV]. The ion with the highest ioniza- 
tion potential (S[IV]) appears to trace the dust continuum best, 
but the lack of more extensive spectral coverage makes it very 
hard to draw any solid conclusions. 
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The various (weak) PAH features have symmetrical inten- 
sity distributions, very similar to each other (Fig. [7]i and con- 
sistent with the 8.0 //m image which shows the stronger emis- 
sion in southwest of the X-ray contours. Fig. [7] also shows the 
[NeII]/[SIV] ratio. The excitation potentials of Nell and SIV are 
21.6 and 34.8 eV respectively. The ratio of the corresponding 
ionic lines [NeII]/[SIV] can be used as a tracer of the hardness 
of the interstellar radiation field (ISRF) in a similar way as the 
ratio [NeII]/[NIII] (see for exemple Gi veon et al.|2002] >: a lower 
ratio corresponds to a harder ISRF. If PAHs were destroyed by 
FUV photons (e.g. |Madden et al.|[2 006), PAH intensities should 
be correlated with the [NeII]/[SIV] ratio. 

Inspection of Fig. [7] reveals a clear trend: the PAH profiles 
peak where the [NeII]/[SIV] ratio dives (+20" W) probably indi- 
cating the edge of an ionized region where PAHs are destroyed. 
We refer for a more extensive and detailed discussion of PAHs 
in the whole 30 Doradus region to |Bernard-S alas| ( [2008) > . 

4.3. Conditions in the extended dust cloud 

Fig. [5] depicts the spatial intensity variation of spectral features 
along the significantly longer SE-NW LL2 slit. The slit extends 
from a region with weak IR emission and no X-rays to the north- 
west, fully crossing the Ha and X-ray emitting SNR (Fig. [2] 
left-hand panel). The compact source J05375027-691 1071 is 
brighter than its surroundings in the continuum, and in the [SIII] 
and [Fell] lines, but not in the [Nelll] line. Beyond this source, 
the distributions of [SIII], [Nell] and the continuum are simi- 
lar and peak about 50 arcsec (12 pc projected distance) to its 
northwest. The [Fell] from the diffuse cloud is relatively sharply 
peaked. Although uncertain, this peak may indicate the presence 
of shocked gas at an SNR/dust cloud interface. As the Ne[III] 
and [SIII] intensity distributions peak closest to the early O stars 
in LH 99 (see Fig.|2]and Sec. 4.3), they seem to reflect the radia- 
tive effects of these. 

The ISO SWS and LWS instruments have provided spec- 
troscopy from differently-sized regions (Vermeij et al. 2002) 
roughly covering the 03 stars marked in Fig.|2| The LL2 slit 
misses the SWS aperture completely, but does overlap with part 
of the LWS aperture. The [Nelll] 15.6/36.0 /mi, [SIII] 18.7/33.5 
//m, and [OIII] 51.8/88.4 pm ratios (6.73, 0.40, 0.56 resp.) are 
all in the low-density limit, indicating n e < 100 cm 4 . The ISO 
spectra also show [Fell] 26.0 //m and [Sill] 34.8 /mi line emis- 
sion but not sufficiently dominant to be ascribed unequivocally 
to shock excitation. More importantly, the ISO data by Vermeij 
etaLl ( |2002] l imply ratios [SIV]/[SIII] = 0.25 and [NeIII]/[NeII] 
= 1 .74. These values are also in excellent agreement with the re- 
lation established by Martin-Hernandez et al. (2002), and place 
the cloud somewhat closer to the bright LMC nebulae. Although 



there are complications briefly discussed by Martin-Hernandez 
|et al. ( |2002| l, the ratios, taken at face value, suggest excitation by 
05-06 stars rather than by 03 stars. 

The spectral energy distribution of the extended infrared 
cloud is shown in Fig. [8] The emission peak is consistent with 
a modified black-body of temperature T = 29 K (B(v, T) x v' 3 , 
with /3 = 2). Integration of the emission yields a flux of 1.1 - 
1.5 x 10~ u W irT 2 implying a luminosity L = 0.8-1.1 x 10 6 
Lq. There is a second hot dust component, with a tempera- 
ture T ~ 230 K, and a luminosity (excluding emission from 
J05375027-6911071) L = 4 - 5 x 10 5 L Q . In the 12 fim IRAS 
band the compact source has a strong silicate absorption at 10 
fim, which makes not possible the separation of the compact 
source and extended cloud contributions. Before fitting the hot 
dust component we thus removed the 12 fim IRAS point, for 



which the flux density is in fact 30 times lower than for the 
neighboring bands. Because of the very strong temperature de- 
pendence of dust emissivity, the mass in the hot dust component 
represents only a minute fraction of the total mass. The three 03 
stars dominating the LH 99 population are projected onto the 
northeastern edge of the infrared cloud as depicted at 24 fim in 
Fig.[T] To their southwest, i.e. closer to the brightest part of the 
nebula, there are at least 2 05 stars, 5 06 stars, and 5 07 stars 
(cf. Schild ~& Testor|1992| ). Using the luminosity scale by |Vacca| 
|et al.| ( |1996) , we find that these stars add another L(O5-07) = 
4.1 X 10" Lq to the luminosity L(03) = 3.6 x 10 6 L Q already pro- 
vided by the brightest stars. The hardness of the resulting mean 
radiation field should therefore not be very different from the one 
estimated above. 

Emission from the radio source appears to avoid the dark 
cloud (see Fig. 4 in Dic kel et al.| |T994). This is not expected if 
the dust cloud is only a foreground object, and it therefore in- 
dicates physical contact between the radio (SNR) and infrared 
(molecular cloud) sources respectively. However, the infrared 
cloud reradiates at most 15% of the total photon output of the 
association members, requiring it to cover only about a steradian 
as seen from the average O star. LH 99 should thus be deemed 
to be quite capable of providing the energy to power the infrared 
cloud, and there seems to be no need to invoke energy inputs 
provided by an SNR impact. 

|Chu et ak] ( |1992| l have suggested that no more than about 
20% of the radio emission from MC 69 is thermal in origin, so 



that S 5GHz(th) < 0.45 Jy (cf.|Lazendic et al. 2000|). This requires 



a maximum Lyman continuum photon flux Nl < 1 .0 x 10 50 s 1 . 
The three somewhat excentric 03 stars provide a combined flux 
N L = 2.5 x 10 50 s _1 , and the more embedded later-type O stars 
an almost equal N h = 2.4 x 10 50 s _1 . Thus, the UV photon flux 
likewise should be sufficient to provide for the ionization of the 
HII nebula associated with the SNR and dust cloud. 



4.4. The SNR revisited 

The actual location of the cloud with respect to the remnant is 
still poorly determined. The data presented here have revealed 
no evidence for shocked material, although the dense cloud of 
dust and molecules does appear to be quite close to the SNR. Is 
the absence of such impact (shock) signatures consistent with the 
dynamical evolution of the remnant? We assume that the rem- 
nant expansion is in the Sedov phase (Wang & Gotthelf 1998), 
so that the temporal evolution of the spherical blast wave with 
radius R$ , produced by the supernova explosion into a medium 
of uniform density po, is described by (Sedov 1959; McCray & 
|Snow|1979): 



R s = 1.15(£ /po) 1/3 f 



1/5 fi/S 



13 (E 51 /n ) l ' 3 t: 



1/5 J/5 



pc. 



We assume an initial energy E st _ Eg/(10 51 erg s ') = 1, a rem- 



nant age f 4 = f/(10 4 yr) = 0.5 (Marsha ll et al.|l998|, and a rem 
nant radius R s 1 5 pc (|Chu et al.||1992[ |Dickel et al.p9"94 



Lazendic et al. 2000 Chen et al.|20 06 ). Thus, the present size of 
the supernova remnant requires an initial density of no more than 
«o = 0.12 cirT 3 . This rules out the possibility of N157B having 
expanded into an interstellar gas cloud of any significant density, 
and confirms the notion that expansion occurred inside a wind- 



blown cavity cleared out by the supernova progenitor ( Chu et al 



19921 ) In the ir X-ray study of 30 Doradus and N157B, |Townsley 
et al.| ( |2006] ) present a composite X-ray /Ha/mid-IR image (their 
Fig. 14) that indeed shows a shell of ionized gas and warm dust 
surrounding the SNR. To verify whether a suitable low-density 
cavity could have been generated by the supernova progenitor, 
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we calculate its expected size by fMcCray & Snow 1979| and 
references therein): 

1/5 -1/5.3/5 

~ 27 L 3 ' 6 n t 6 ' pc 
where R$ is the cavity radius, the constant wind luminos- 
ity in units of 10 36 erg s _1 , «o the ambient gas particle den- 
sity in cirT 3 , and % the time that the wind has been blow- 
ing in millions of years. The wind luminosity is given by: 
L 36 * 0.3(M„/20M G ) 2 - 3 erg s _I We adopt a M t = 25 M (cor- 
responding to an 08-09 star), intermediate between the values 
pre dicted by the |Thielemann et al.| ( fl996] l, |Heger et al] ( |2003] l 
and |Woosley & Weaver ( 1995| l models so that = 0.5. The 
relevant time t(, = 3.2 is the enti re lifetime of the progenitor, 



derived from e.g.|Karttunen|(|1995|>: t ~ (MJM @ T 2 - 5 x 10 10 yr. 



The observed radius R$ of about 15 pc implies an average den- 
sity of «o = 310 cirT 3 for the ambient medium into which the 
windblown shell has expanded. Such a value is compatible with 
the densities of a few hundred per cc that we estimated in the 
previous section for the clouds in the region. Even the larger ra- 
dius Rs ~ 22 pc gleaned from |Townsley et al. ( |2006 1 can be 
accommodated as it still yields «o = 50 cm" 3 . This situation 
changes rapidly, however, if we consider a more massive super- 
nova progenitor. For instance, for M, = 50 Mo and = 15 pc 
we have L36 = 2.5 and fe = 0.6, and find an average ambient 
density of only no = 9 cirT 3 . 

Thus, we conclude that the lack of clear indications for an 
SNR impact on the dust cloud seen in the IRAC and MIPS im- 
ages is consistent with the supernova explosion of a moderately 
massive star. Only now would the SNR expanding in the cav- 
ity begin to overtake the windblown shell produced by the star 
over its lifetime. The apparent lack of shocks is not easily recon- 
cilable with a supernova progenitor as massive as the presently 
most luminous members of the LH 99 association. This confirms 
an independent conclusion by Chen et al. (2006) that the progen- 
itor should have been in the narrow mass range M = 20-25 
Mo. 



5. Conclusions 

From an analysis of Spitzer photometric mapping and spec- 
troscopy of the SNR-dominated region N157B in the LMC, we 
find that: 

- There is no evidence of an infrared counterpart to the super- 
nova remnant in the IRAC and MIPS images. 

- The infrared emission is dominated by a cloud of dust and 
molecular gas adjacent to the remnant, containing the com- 
pact 2MASS source J05375027-691 1071. 

- The object J05375027-691 1071 has a diameter of about 3 pc, 
an electron-density of 100-250 cm 4 , and is photo-ionized by 
an 08-09 star. It is probably an open HII blister structure, 
seen from the back. 

- In spite of the projected overlap between the SNR X-ray 
emission and the infrared cloud, there is at best very marginal 
evidence of shocked gas, while almost all data suggest 
photo-ionization and photon-heating to be the mechanisms 
dominating the infrared cloud. 

- The extended dust cloud is associated with ionized emission 
of a density of typically 100 cm" 3 , presumably at the edges 
of a denser molecular cloud. As the extended dust reradiates 
only about 10 per cent of the luminosity of the 15 brightest 
and nearby O stars in the LH 99 OB association, these stars 
are sufficient to explain the heating of the dust cloud. 

- The absence of clear evidence of shocks implies that at 
present the molecular/dust cloud is not significantly im- 



pacted by the remnant. This suggests that the supernova pro- 
genitor was a moderately massive star of mass M » 25 M . 
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